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ABSTRACT 


The so-called drag-based model (DBM) simulates analytically the propaga¬ 
tion of coronal mass ejections (CMEs) in interplanetary space and allows the 
prediction of their arrival times and impact speeds at any point in the heliosphere 
(“target”). The DBM is based on the assumption that beyond a distance of about 
20 solar radii from the Sun, the dominant force acting on CMEs is the “aerody¬ 
namic” drag force. In the standard form of DBM, the user provisionally chooses 
values for the model input parameters, by which the kinematics of the CME over 
the entire Sun-“target” distance range is defined. The choice of model input pa¬ 
rameters is usually based on several previously undertaken statistical studies. In 
other words, the model is used by ad hoc implementation of statistics-based val¬ 
ues of the input parameters, which are not necessarily appropriate for the CME 
under study. Furthermore, such a procedure lacks quantitative information on 
how well the simulation reproduces the coronagraphically observed kinematics of 
the CME, and thus does not provide an estimate of the reliability of the arrival 
prediction. In this paper we advance the DBM by adopting it in a form that 
employs the CME observations over a given distance range to evaluate the most 
suitable model input parameters for a given CME by means of the least-squares 
fitting. Furthermore, the new version of the model automatically responds to any 
significant change of the conditions in the ambient medium (solar wind speed, 
density, CME-CME interactions, etc.) by changing the model input parameters 
according to changes in the CME kinematics. The advanced DBM is shaped 
in a form that can be readily employed in an operational system for real-time 
space-weather forecasting by promptly adjusting to a successively expanding ob¬ 
servational dataset, thus providing a successively improving prediction of the 
CME arrival. 
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Subject headings: Sun: corona — Sun: coronal mass ejections — solar wind: 
solar-terrestrial relations — magnetohydrodynamics — methods: analytical - 
methods: numerical 
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Introduction 


Eruptive processes in the solar atmosphere, particularly coronal mass ejections 
(CMEs), strongly influence the physical state of the heliosphere and the terrestrial space 
environment. CMEs represent eruptive restructuring of the global coronal magnetic field, 
where the eruption itself is caused by a loss of equilibrium of the pre-eruptive magnetic field 
structure. The stability of the structure depends on the amount of energy stored in the 
magnetic field, whereas the CME itself is driven by the Lorentz force. The dynamics of the 
instability depends on the magnetic-flux conservation and inductive effects, which cause the 
cessation of the Lorentz force. Eventually, the magnetohydrodynamic (MHD) drag becomes 
a dominant factor in the CME dynamics. The drag is a consequence of collisionless transfer 
of momentum and energy between the CME and the ambient solar wind by MHD waves 


(Cargill 2004). 


In the present paper, we develop a method that provides the observations-driven 
adjustment of the input parameters of the so-called drag-based model (hereinafter, DBM), 
which describes the CME propagation in the interplanetary space by considering the “drag” 
force (for details see Vrsnak et al. (2013), and references therein). The “drag” force depends 
on the relative speed of the ejection and the solar wind; in a collisionless environment the 
acceleration can be expressed as a = —7 (u — w) \v — w\, where 7 is the “drag parameter”, a 
and v refer to the instantaneous acceleration and speed of the ejection, whereas w represents 
the ambient solar wind speed (Vrsnak 2001, Cargill 2004| Owens & Cargill 2004, Vrsnak 


& Zic 

2007, 

Borgazzi et al. 

2009; 

Lara & Borgazzi 

2009, 

Vrsnak et al. 

2010 , 

Vrsnak et al. 


2013). Furthermore, the previously used DBM with constant 7 and w parameters (Vrsnak 


et al. 2013) is extended into a more general form, allowing variable 7 (r) and w(r). In 


the DBM the CME is represented by the cone shape, where each element of the CME’s 
leading edge is defined by its position relative to the CME tip. The parameters 7 and w 
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represent the most sensitive elements of the DBM and play the main role in the drag-based 
simulation of heliospheric CME propagation. Consequently, their evaluation represents the 
central issue in DBM-based space-weather forecasting. 


The paper is focused on the theoretical elaboration of finding values of the DBM 
parameters that give the smallest difference between the DBM-based kinematics and the 
CME kinematics as derived from observational data. The observational measurements could 
be derived from coronagraphic and heliospheric imaging data using several methods based 


on certain assumptions (e.g., fixed (p for small CMEs: see Sheeley et al. 1999; Rouillard 


et al. 2008 or harmonic mean for large CMEs: see Lugaz et al. 2009). The presented fitting 


method opens the possibility of an “automatic” evaluation of the most appropriate DBM 
input parameters from observational data available for a particular event. The application 
and validation of the proposed method will be presented in a follow-up paper employing 
detailed coronal and heliospheric observations of one slow and one fast CME. 


2. General description of the drag-based model 


2.1. The drag force 


In interplanetary space the CME motion is governed by the Lorentz force F f, gravity 


F g , and the MHD analog of the aerodynamic drag F ( \ (Vrsnak 2006). The net CME force 
can be expressed as: 

F = F l - F g + F d . (1) 


At heliocentric distances beyond R > 15, the MHD drag becomes a dominant force (Vrsnak 


et al. 2009]), so the CME motion is basically influenced solely by the Fa term of the force 
Equation (jlj. 


Generally, the “drag” interaction between the solar wind and the CME in interplanetary 
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space can be described in various ways. In this paper we consider the “drag” force of the 
form: 

Fd = — Cd Ap (v — w) |u — w \. ( 2 ) 

where Cd refers to the dimensionless drag coefficient, A is the cross section of the CME, p 
represents the ambient solar wind density, and (v — w ) is the velocity difference between 


the CME and the solar wind (Chen 1989, Chen & Garren 1993, Cargill et al. 1996, 2000 


Cargill & Schmidt 2002 Vrsnak & Gopalswamy 2002 Vrsnak et al. 2004 2009, 2010[ ). 


Following the numerical MHD simulations by Cargill (2004), and under the assumption 


that the CME structure does not change, we expect that the drag coefficient Cd varies 
slowly with radial distance and is approximately equal to 1 for the heliocentric distances 
beyond 15 solar radii, particularly in the case of dense CMEs. The mass density of a CME 
lies in the range of 12.7 — 13.5 g/ tq, with the most dense events occurring during the solar 
maximum (see Vourlidas et al. (2011)). In this respect, we define dense CMEs as events 


with a mass density exceeding 13.2g/r|,. 


The CME acceleration, caused by the MHD “drag” (Vrsnak et al. 2009), can be written 
in a simple form using Equation (|2]) : 


ad = —7 {y — w) \v — w\, 
where the parameter 7 is defined by 


( 3 ) 


7 = Cd 


Ap 

M' 


( 4 ) 


The parameter 7 is inversely proportional to the total CME mass M, which consists of the 
initial mass and the so-called virtual mass that piles up as the CME expands in the inner 
heliosphere. Observations indicate that beyond heliocentric distances of several solar radii 


the total mass becomes approximately constant (Bein et al. 2013), implying that the mass 


pile-up becomes balanced by the mass loss (Vrsnak & Zic 2007, Vrsnak et al. 2013). 
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2.2. Ambient density and solar wind speed 


For the ambient density po( r ) = ^p '^o (?"), where no(r) is unperturbed particle density 


and m p is proton mass, the empirical no(r) model proposed by Leblanc et al. (1998) is 
applied (referred to in the following as LDB, after Leblanc, Dulk, Bougeret). The CME 
cross section A(r) depends on the geometrical shape of the CME; hereafter the cone 


representation is employed (Fisher & Munro |1984 ; Xie et ah 2004; Schwenn et ah 2005). 
The parameter 7 (r) in Equation ^ defines the effectiveness of the drag force and depends 
on both the CME and the ambient solar wind. 


At distances beyond R> 15 (R = r/r Q , where r Q is the solar radius), the terms oc 
and oc R ~ 6 in the LDB expression for no(R) can be neglected. However, for the purposes of 
completeness and model development, in this paper the complete LDB density expression 
is applied: 


, x-,\ k 2 kd kp. 

na(R) = W + W + 


( 5 ) 


which is valid for radial distances R > 1.8. Coefficients k 2 , k 4 and k$ read k 2 = 3.3 x 10 5 cm 


-3 


fc 4 = 4.1 x 10 6 cm 3 and k 6 = 8.0 x 10 7 cm 3 (Leblanc et al. 1998). 


The background solar wind is taken to be approximately stationary and isotropic, so 
from flux conservation, drio/dt + V ■ (now 0 ) = 0, the solar wind speed must satisfy the 
expression: 

/m f, , h/k 2 h/k 2 \~ l 

w 0 (R) = ( 1 + -^ 2 “ + -jp- I , (6) 

where is the asymptotic solar wind speed, i.e., w^ = lim^oo w 0 (R) = const. At small 
heliocentric distances the solar wind rarefies at a rate larger than R~ 2 , so the wind speed 
has to rise according to the continuity equation. Figure [l] shows the radial dependences 
of the normalized drag parameter 7 , solar wind speed Wq and density Hq. The presented 




















ratios r ){R)/ r ) 00 and w 0 (R)/w oo are normalized by asymptotic values ( 7 ^ and w^), and the 
ratio n 0 (R) /uiau by the density value Uiau at 1 AU. Note also that the value of 7 close to 
the Sun is for an order of magnitude larger than it is at large distances. As can be seen 
from Figure [lj one finds that the dependences Wq(R), as well as 'y(R), become practically 
constant beyond R > 15. Thus, the asymptotic values of w 0 and 7 are approximately equal 
to the values at 1 AU, i.e., ~ wiau and 7 oo ~ 7 iau- A similar simplihcation was used 


in previous papers (Vrsnak & Zic 2007, Vrsnak et al 


2010 


Vrsnak et al. 2013) where the 


unperturbed solar wind speed and parameter -7 functions had constant values, w 0 (R) = w^ 
and 7 (R) = 700 , for all radial distances. 

The solar wind speed Wq(R) can be additionally modified by including a specific 
perturbation w p (R ) on top of the described undisturbed background to reproduce a 
particular situation in a given event. For example, in some cases the CME travels in 
interplanetary space through a region of locally enhanced or decreased solar-wind density 


(Temmer et al. 2011, 2012, Maricic et al. 2014, Rollett et al. 2014). In such a case the 
additional w p (R) term should describe the associated solar wind speed perturbation in 
the region between the heliocentric distances R\ and R 2 - Under these assumptions, the 
perturbed solar wind speed is defined as 


w{R) = 


w 0 (R ) + w p (R), Ri < R < R 2 
wq(R), otherwise 

where Wo(R) represents the unperturbed solar wind speed (see Equation 
perturbed density induced by the solar wind term w p in w(R ), 


(7) 


The 


n{R) 


h '2 ^00 

R?w(R)’ 


( 8 ) 


follows from flux conservation, i.e., n(R) = lim Roo ^. 00 [?r 0 (i? 00 ) R 2 ^ Wq(R 00 )/w(R) R 2 ]. 
The perturbation is assumed to be localized over a finite region (i.e., inside the interval 
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Fig. 1.— Radial dependence of the solar-wind density no normalized by the density niAu 
at 1 AU (n 0 /niAu; solid line), shown together with the solar wind speed wq and the drag 
parameter ratio 7 normalized by their asymptotic values and (dashed and dotted 
lines, respectively). 



Ri < R < R ‘)), whereas the unperturbed expressions for the density, Equation ([ 5 ]), and 
solar wind speed, Equation ([6]), are valid otherwise. 

The solar wind perturbation above is described by defining the wind speed w p (R) from 
which follows the density prohle n(R); however, if the case study requires, the perturbation 
could be performed in the opposite way, firstly defining the density perturbation n p (R) and 
afterward evaluating the solar wind speed expression w(R). 


1AU 
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2.3. Drag parameter 7 

In the case of the CME propagation where the “drag” force is dominant, the equation 
of motion, Equation (J3]) , transforms to 


m 


- 7 (R) R(t ) - m(R) R(t) - w(R) 


( 9 ) 


The expression for 7 given by Equation Q shows that the drag is more effective if the 
ambient density is high, if the CME is light, and if the CME cross section is large. I 11 the 
present version of the DBM the effective CME cross section is defined by employing the 


CME “cone model” (Fisher & Munro 1984, Xie et ah 2004; Schwenn et al. 2005). In this 
presentation, the cross-sectional area is given by A = R 2 tan 2 A/(l + tan A) 2 , where 
A is the CME half-width (see the Appendix). Note that CMEs could be represented by 
a variety of geometrical representations; for examples of commonly used geometries see 


Schwenn et al. (2005), Thernisien et al. (2006), Lugaz et al. (2010), Thernisien (2011), 


Davies et al. (2012), and the references therein. Note that in the previous form of DBM 


with constant w and 7 (Vrsnak & Zic 2007, Vrsnak et al. 


2010 


Vrsnak et al. 2013) the area 


was approximated by A & nr^R 2 A 2 . However, this did not have a direct influence on 
the calculated CME kinematics, since the value of A was already incorporated within the 
presumed parameter 7 . 


Taking into account the definition of solar wind speed, Equation (JT]) , 7 includes even 
the cases of perturbed solar wind, i.e., when the term w p (R) is taken into account: 


7 (R) = 7c 




w{R) 


( 10 ) 


Since the asymptotic value of the solar wind speed at large heliocentric distances 
(R —>■ 00 ) is Woo (see Figure [I]), evidently 7 (R) asymptotically acquires value 700 likewise, 
i.e., 7 oo = lim^oo 7 ( 7 ?). In space-weather forecasting it became a practice to use a 
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dimensionless variant T, defined by 7^ = T x 10 7 km 1 . 

An interesting consequence follows from the dependence of the parameter T on the 
CME’s geometrical shape and its properties. In the case when the effective cross section 
A is proportional to the distance squared, A oc R 2 , T could be generally calculated using 
T = The expression /(A) depends on the CME geometrical shape of the CME used 

in the model, which in our case is given by /(A) = nm p r q [tan A/ (1 + tan A)] 2 x 10' km = 
8.4 x 10 12 [tan A/ (1 + tanA)] 2 kg (see the Appendix). If the observations provide the 
CME half-width angle A, and observed kinematics provide the value of T, the CME mass 
M can be roughly estimated by employing M = 8.4 x 10 12 [tan A/ (1 + tanA)] 2 r _1 kg, 
where A is expressed in radians and M in kg. The same holds for the opposite situation: 
knowing the mass M one can estimate the angular half-width A from the value of T. The 
presumed geometrical shape of the CME affects the estimation of cross-sectional area, and 
consequently is important in evaluation of the unknown properties (M or A) of the CME. 
Thus, the best suited choice of the geometrical model to the CME observational properties 
improves the accuracy of the CME mass or half-width evaluation. 


3. Model/observations fitting 

In the following, the procedure of finding the values of any unknown DBM parameters 
is described. The drag parameter T, the background solar wind speed w^, and the modified 
initial CME radial distance Rq and speed vq are adjusted iteratively by minimizing the 
deviation of the model kinematics from the observed one. The process sequentially alters 
the DBM parameters in order to minimize the quadratic deviation (the sum of squared 
“errors” or residuals) between observational and DBM-calculated speeds: 

N 

E(T, w oo; R 0 , t’o) = ^ [vj ~ v({r, wR 0 , u 0 }, Ri)] 2 ■ 

i= 0 


( 11 ) 
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The observational distance-speed values are written as Ri and v i: while the adjusted 
kinematic curve v(R), dependent on the model input parameters T, w 00 and the initial-state 
parameters Rq, vq, is designated as u({T, w^, i? 0 , uo}> R)- The initial-state parameters 
depend on the presumed geometrical representation of a CME; therefore for the geometrical 
option presented in the Appendix, the initial-state parameters Rq and y 0 of Equation ( JTT| 
can have the CME tip values Ro(t 0 ), i’o(to), or the flank values R v (t o), v v (to) at the 
initial time to, depending on the observer’s location (see the Appendix). Unknown input 
parameters are found by successively solving the equation of motion, Equation ([9]), within 
the parameter domain. In practice we use a more appropriate form of Equation (|9|, which 
reads 


v(R) 


d v(R) 
d R 


= - 7 (R) [v(R) - w(R)] \v(R) - w(R )\. 


( 12 ) 


The variation of the DBM parameters seeks the minimal value A mill of Equation 0 - 
The presented method is basically a modified successive multiparametric variation that 


includes solving of the differential equation of motion, Equation (12), and least-squares 
fitting (hereafter, LSF) to the observational {(R 0 , v 0 ), ■ ■ ■, (Rn, vn)} dataset. Different 
approaches could be used in the numerical fitting. For example, the computation could be 
performed by starting with arbitrary DBM values based on which optimal values are found, 
or by numerically seeking the minimum of Equation ( [TT| ) within a physically meaningful 
DBM-parameter domain (Motulsky & Ransnas 1987). The meaningful parameter-domain 
restriction could be also included in the firstly mentioned approach to speed up the 
process of finding the E m[n . In the end, the minimal quadratic deviation gives the best 
input-parameter set {T, w^; Rq, Uo} for the specific observational event. Furthermore, 
kinematic curves, such as a(R), a(t), v(R), v(t), and R(t) are automatically available from 
the calculated parameters. Consequently, this directly provides the CME transit time r, 
defined as the time the CME takes to arrive at a prescribed location, as well as the “impact” 
velocity v T . 
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It is instructive to employ statistical analysis and to express the “goodness” of fit in 
the form of several statistical quantities. The first is the standard deviation (or the rms), 
which represents the average deviation between observed v t and calculated v(Ri) data: 



and gives data dispersion in velocity units (i.e. kms : ). Notice that the (TV + 1) is the total 
number of observational datapoint samples {(i?o, ^o)> •••■., (Rn,vn)}- 

The “goodness” could be graphically presented in the form of a residual plot. A 
residual plot shows the differences (or residuals) between each measured v t value and the 
value calculated from the estimated curve v(Ri), i.e., — v(Ri)]. The residuals should not 

have a systematic dependence on R values (the abscissa values) and should have a random 
scattering. Any clustering of residuals in the plot indicates that the prediction curve follows 
a systematic-error pattern and that the fit is not appropriate. 

The next relevant criterion of scattering between observed values Vi(Ri) and calculated 
v(R) is the coefficient of variation, which is defined by 

c v = -- 100%, (14) 

v 

where v = ^2^ =0 v(Ri) / (N + 1) is the mean of calculated values from the calculated set 
{»’«)}• 

Lastly, the coefficient of determination is defined by 

R 2 = 1 - ^=°j Vi ~ v ( Ri ^ 7 0 < R} < 1 (15) 

E t=0 k - v \ 

The “quality” of the estimated DBM parameters increases, i.e., the calculated kinematic 
curve fits observational data better, as a and c v decrease and reach minimal values cr m i n and 
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Cvmim respectively. On the other hand, 1Z 2 becomes close to 1 as fit gets better (Motnlsky 
& Ransnas 1987:). 


For demonstration purposes we have chosen the observational dataset of Event 1 


described in Temmer et al. (2011) and we have applied the LSF-DBM method to 
simulate the CME propagation. The event started on 2008 June 1 at ~ 21 UT and the 
propagation data {Ri,i>i} with associated errors are derived from STEREO coronagraphic 


and heliospheric image data using the constrained harmonic mean method (see Rollett 


et al. 2012). For more details we refer to Temmer et al. (2011). The final result of the DBM 


fitting is presented in Figure [2j where panel |2^a) presents the velocity-distance profile of the 
calculated CME kinematics (blue solid line and accompanying shaded error area), as well as 
the estimated solar wind speed (green dashed line), together with the observational dataset 
(black circles) and its error bars. In the bottom panel Qb)) we present the residulas, i.e. ; 
the relative difference pty — v(Ri)]/vi between observational, Vi, and calculated, v(Ri), CME 
velocities, relative to the velocity v t . The standard deviation of the observed dataset is 
<t 0 = \J^2iLo /(-^ + 1) = 42.81 kms -1 , where e* represents the half-error bar value for 
each measurement of velocity ly. The LSF-DBM technique produced the fit with the DBM 
parameters T = 2.84, = 433.04kms -1 , v 0 = 229.50kms -1 , R 0 = 14.17r Q , accompanied 

by the minimal standard deviation and the coefficients of variation and determination, 
u min = 29.87 km s -1 , c V min = 7.50%, 77 2 = 0.67, respectively. In Figure |2](a) the “average” 
errors of the fit, spanning values frim v(R) — <r min to v(R) + cr min , are drawn as the blue 
shaded area in the vicinity of the DBM kinematic curve v(R). Evidently, the fitted standard 
deviation a m \ n is much smaller than the observed <r 0 , showing that the LSF-DBM produced 
a satisfactory fit within the range of the “average observational error” cr 0 . 
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Fig. 2.— Application of the LSF-DBM method on the dataset extracted from Event 1 of the 


paper by Tennner et al. (2011). (a) Radial dependence of the LSF-DBM estimated kinematic 
curve v(R ) (solid blue line with error shown as the blue shaded area) on the observed CME 
speed values (black circles, with error bars), and numerically calculated solar wind speed w 
(dashed green) based on parameter fitting, (b) Radial dependance of residuals, Vi — v(Ri ), 
between observational Vi and DBM-calculated v(Ri) speeds of CME (relative to observed 
CME speed U;). 
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4. Discussion and conclusion 


We have presented an extension of the DBM that is intended for use in automatic 
forecasting of CME arrival and impact at an arbitrary heliospheric position. The extension 
consists of optimizing the DBM input parameters based on the sequential variation and 
determination of the minimal standard deviation (cr min ), the minimal coefficient of variation 
(Cy min )i and the coefficient of determination ( 1Z 2 ) from observational data. The mentioned 
statistical quantities represent an estimate of “goodness” of the DBM fit to observational 
data and consequently define the reliability of the arrival time and impact speed prediction. 
The presented LSF-DBM modification opens an opportunity for implementation in 
real-time space-weather forecasting tools and alerting systems for CME impacts on Earth 
(or any heliospheric “target” of interest). The novel approach is based on real-time 
data-driven DBM-parameter optimization that iteratively improves the accuracy of CME 
kinematics in the heliosphere. 


The accuracy of the real-time forecasting increases as the observational dataset 
successively becomes larger, i.e., as the CME is tracked to larger heliospheric distances 


(Davis et al. 2010). For example, we can imagine a hypothetical case study of a CME 
launched from a region close to the solar disc center at a specific time t — 0. In this case the 
CME is directed toward Earth and we can extract information about a current in situ solar 
wind speed at the Earth. On the other hand, the CME is traced in real time during 
its propagation throughout the heliosphere, so by using the most likely CME geometry 
the observational data can be transformed to get the distance-speed ( Ri,Vi ) data. Every 
time when a new (i?j,Uj) datapoint become available, the LSF procedure estimates a new 
set of DBM parameters required for updating the DBM forecast of CME arrival. As the 
dataset expands, our “impact prediction” becomes more reliable. In this respect it should be 
noted that an L5 mission is urgently needed to advance the performance of such forecasting 
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methods. We note that this example is a quite simplified case in which the LSF-DBM 
method could be used. 

There are several drawbacks of the described procedure, e.g., the estimation of the 
model input parameter and the related forecasting are highly dependent on the quality of 
the observational dataset. The input required for the DBM fitting procedure is the observed 
set of values for speed and distance of the CME frontal part as observed along the ecliptic 
plane. Several methods exist to derive those quantities, and we just mention briefly some 
possibilities. The propagation direction might be simply estimated from the CME associated 
source region, assuming radial propagation. In fact, knowing the propagation direction 
would enable one to derive the 3D CME kinematics from single spacecraft observations, 
such as, e.g., from STEREO heliospheric image data assuming a certain CME width (e.g., 
fixed </> for small CMEs: see Sheeley et al. 1999} Rouillard et al.|2008 or harmonic mean 


for large CMEs: see Lugaz et al. 2009). Using stereoscopic data, triangulation methods 


could be used that also provide the required input (coronagraphic held of view: see Mierla 


et al. 2010, interplanetary space: see Liu et al. 2010). The uncertainty of measurements 


is automatically forwarded to the estimated model parameters, and consequently, to the 
arrival time prediction. 

Another serious drawback lies in the fact that the employed observational data include 


the distance range where the CME is still driven by the Lorentz force (Gallagher et al 


2003). In such a situation, the DBM fails in its fundamental concept because the Lorentz 


force is excluded from the modeling, i.e., only the drag force governs the CME propagation. 
However, it should be noted that even in such a case, the DBM kinematical curve might 
fit the observational data nicely due to the fact that the statistical weight dominantly 
comes from larger heliospheric distances, where the Lorentz force should be negligible. For 
example, if the observational dataset used in the modeling consists of only a few low-height 
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measurements and a more abundant subset of measurements at larger heliocentric distances, 
the larger drag-dominated dataset “overweights” the smaller Lorentz-driven subset, so the 
latter effect becomes negligible. 


The LSF and the DBM could be used in an opposite way to that previously discussed, 
for example to estimate the solar wind speed w^ at large heliospheric distances (R > 15 r Q ), 
in the regions where in situ measurements are not available. Moreover, measuring ( Ri,Vi ) 
LSF straightforwardly gives the solar wind speed, w^. Using Equations (Jr]) and ([8]) we 
can then roughly calculate w(R) and n(R) for any heliocentric distance, R. Additionally, 
as the LSF estimates the complete set of DBM parameters, {T, Ro, no}, i n situations 
when the measurements are not very confident and have a high uncertainty, we could apply 
the LSF method and correct, for example, the low-coronal initial position and the velocity 
of a CME. However, the unknown DBM parameters are more reliably estimated as more 
parameters are directly given from the observations, and if stereoscopic observations are 
conducted in an appropriate manner to provide reliable deprojected (R i} Vi) values. 


The LSF-DBM could be further applied in a case when a CME meets various 
heliospheric “obstacles” during its propagation. The probability of an interaction between 
two consecutive CMEs is very high in the heliosphere, since on average several CMEs are 


observed per day with different kinematics and velocities (St. Cyr et al. 2000, Gopalswamy 


2006). The interaction takes place when the later and faster CME catches an earlier 


and slower one (Temmer et al. 2012; Maricic et al. 2014). By inspecting the CME’s 
behavior and surrounding ambient conditions, the LSF-DBM procedure could be used for 
a “segmented-distance” application. For example, the CME trajectory could be divided 
into several parts dependent on the CME behavior, i.e., divided into regions before the 
CME-CME interaction and the region after the interaction. In that way the forecasting 
of the CME arrival at a given “target” could be acquired by applying sequentially the 
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LSF-DBM technique on each trajectory interval (e.g., Temmer et ah 2011, 2012, Maricic 
et al. 2014; Rollett et al. 2014]). 


Actually, numerical computation requires arbitrary initial DBM-parameter entries 
around which the LSF procedure searches for the best result. Sometimes the problem arises 
when, in the proximity of starting DBM entries, the numerical LSF finds multiple cr min 
minima inside the parameter domain. The problem could be avoided by carefully studying 


the specific case, or using a different track-fitting method (see, e.g., Mostl & Davies 2013 


who use a constant-velocity approximation) and then reapply the LSF-DBM procedure to 
refine the forecasting. 

The presented generalized DBM is an extension of the model with the assumption 


of a constant 7 (R) and w(R) (Vrsnak et ah 2013), which is not adequate for describing 


low-coronal CME propagation, or kinematics in the spatially perturbed solar wind w(R). 
Finally, the application of the least-squares fitting method coupled with the DBM applied 
to various CME geometries and solar wind models offers an improvement in efficiency and 
accuracy of forecasting CME. 
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A. CME shape used for DBM online tool 

We briefly discuss the general outcome of a DBM calculation and its dependence on 
the presumed shape of the leading edge of the interplanetary CME (ICME). To evaluate 
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the CME cross-sectional area, A, (and therefore h) in this version of the DBM we used 
the geometry presented in Figure [3] (for other frequently used options see, e.g., Sheeley 


|et al.||T999 ; Kahler fe Webb 2007 Lugaz et al. 2009 |Davies et al. 2012 and the references 
therein). In this option the leading edge is considered to be a semicircle, spanning over 
the full angular width of the ICME, 2A. Considering the geometrical relationships between 
various parameters marked in Figure |3j the heliocentric distance R v (t) and the speed v v (t) 
of an element at the angular position p depend on the heliocentric distance of the CME tip, 
Ro(t), the speed of the CME tip v 0 (t), the cone half-width A (which stays constant during 
ICME propagation), and the angle p. Precisely, the relationships between radial distances 
and velocities of the tip and a flank CME element are 


R<p(t) = R 0 (t)F(<p ) 

v v (t) = v 0 (t)F(p ) (Al) 


respectively, where the angular function is the same for both expressions: 


F(f) 


cos p + \J tan 2 A — sin 2 p 
1 + tan A 


(A2) 


The CME expansion is modeled by providing the initial speed Vq and heliocentric distance 
Rq of an arbitrary single point on the CME’s leading edge (e.g., in Figure [3] the leading 
edge segment of the CME tip has the distance -Ro(0) and the speed vo(0)), thus the 
heliocentric distances R v ( 0) and speeds 1^(0) of a certain segment along the leading edge 


with p G [—A, A], follow from Equation (Al). At later time the leading edge evolves 
accordingly, as described in Equation ([9]). 
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Fig. 3.— A cross section of a conical representation of a CME in the ecliptic plane during 
ICME propagation. The simple schematic describes temporal deformation and evolution of 
the CME’s leading edge in time. The initial CME shape at time to is defined from a single 
heliocentric distance measurement R 0 of the CME tip and used in an equation that defines 
the conical geometry. 

In the current versions of the DBM, implemented as a public prognostic online tool at 
http://www.geof.unizg.hr/~tzic/dbm.html, the different options of CME expansions 
are proposed. The prognostic tool forecasts only ICME propagation in the ecliptic plane, 
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as a result of CME initiation at low heliographic latitudes. The geometric setup of the 
latest online DBM version is presented in Figure [3] where the CME frontal part evolves 
in time, i.e., the expansion of the CME’s leading edge is simulated by applying the DBM 
equation of motion, Equation (J9|, on each leading-edge segment independently. The initial 
cross section in the ecliptic plane of the CME shape is constructed by a single /?, o (0) 
measurement of the CME tip element (which lies on the line of the CME’s propagation 
direction and in the ecliptic plane as well) and by the assumption of conical CME geometry 


defined by Equation (Al). The leading edge gradually deforms since different segments are 
immersed at initial time to in different surrounding conditions (described by solar wind 


speed w(R) and y(-R) functions) and have different initial velocities (see Equation (Al)), 
hence the DBM equation of motion results in different radial kinematics. Since the flanks 
move more slowly, and thus in fast ICMEs the drag-deceleration of flanks is weaker whereas 
flank acceleration in slow events is stronger, the variation of speed along the 1CME front 
decreases and the front gradually flattens. Note that such an “independent-element” DBM 
procedure could be equivalently applied to any other presumed initial CME geometry. 
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